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ABSTRACT 


This  paper  considers  a  special  class  of  transportation  problems 
in  which  the  needs  of  each  user  are  to  be  supplied  entirely  by  one  of 
the  available  sources.  We  first  show  that  an  optimum  solution  to  this 
special  transportation  problem  is  a  basic  feasible  solution  to  a  slightly 
different  standard  transportation  problem.  A  branch  and  bound  solution 
procedure  for  finding  the  desired  solution  to  the  latter  la  then  presented 
and  illustrated  with  an  example.  We  then  consider  an  extension  of  this 
problem  by  allowing  the  possibility  of  increasing  (at  a  cost)  the  source 
capacities.  The  problem  formulation  is  shown  to  provide  a  generalization 
to  the  well-known  assignment  problem.  The  solution  procedure  appears 
to  be  relatively  more  efficient  when  the  number  of  uses  greatly  exceeds 


the  number  of  sources. 


1.  INTRODUCTION 


In  a  recent  paper  [4] ,  DeMaio  and  Roveda  consider  a  special  class 
of  transportation  problems  with  a  set  of  sources  I  »  {l ,2, . . . ,i , . . . ,W] 
having  known  capacities  b^  and  a  set  of  uses  J  »  {l,2,...,J,c.,M}  with 
known  demands  r^  for  a  homogeneous  material  (the  and  are  assumed 

to  be  strictly  positive).  The  objective  is  to  minimize  tha  total  traneporta- 
tion  cost  Z  subject  to  the  constraints  that  (i)  each  ucer'a  demand  is 
fulfilled  by  exactly  one  of  the  sources,  and  (ii)  the  total  amount  shipped 
from  each  source  does  not  exceed  its  capacity.  Denoting  by  the  cost 

of  transporting  all  the  r^  units  from  the  ith  source  to  the  jth  uae 
and  defining  to  be  1  or  0  depending  on  whether  or  not  use  j  is  assigned 

to  source  i,  the  problem  is  to 

minimise  Z  *  X  £  c. .x. .  (1) 

lei  jej  3  13 

subject  to  the  constraints; 


T  r.x  <  b.  for  Ul, 

jcJ  J  u  1 

“  x, .  *  1  for  j«J,  and 


k,  ,  *■  0  or  1  for  lei  as  4  (**) 

f  j 

The  authors  of  (4j  present  an  implicit  enumeration  approach  to  aolvtng 
this  problem.  In  Section  2  we  show  that  an  optimal  solution  to  this  problem 
can  be  characterized  as  a  basic  feasible  solution  to  a  eltghtly  modified 
transportation  problem  and  that  such  a  solution  can  be  obtained  by  an  algorithm 
similar  to  the  subtOur  elimination  method  for  solving  traveling  talesman 
problems  t$,  *1.  Since  our  algorithm  utilises  the  underlying  structure  of 
the  transportation  problem,  it  is  believed  to  be  computationally  more  efficient 
tha;>  the  implicit  enuaeration  approach.  As  will  be  seen  in  Section  2  the 
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present  algorithm  appears  to  be  particularly  suitable  whan  the  mother  of 
uses  far  exceeds  the  number  of  sources.  Furthermore,  our  approach  can  be 
easily  extended  to  problems  where  capacity  expansion  for  warehouses  is  a 
possibility.  In  Section  3  we  consider  this  extension  and  provide  other 
practical  applications  covered  by  this  model. 


2.  THE  ALGORITHM 


To  bring  the  problem  (l)-(4)  into  the  standard  transportation  format, 
we  first  make  the  transformations 


y..  =  r  .x.  .  for  id  and  JeJ,  and 
j  ij 

d..  »c,./r.  for  1 el  and  jsJ; 
tj  i  J  J 


(5) 

(6) 


i.e.,  y^  denotes  the  amount  shipped  from  source  i  to  uae  J  at  unit 
cost  d^j.  To  convert  the  inequalities  (2)  into  equations,  we  adopt  the 
usual  procedure  [3]  of  adding  a  slack  use  H  +  1  and  setting 

(1) 

»  0  Cor  id,  and  (8) 

(9) 


l.M+l 


r  «  f  h  «.  «*■ 

M+L  i  “ 

id  jSj 


J' 


The  problem  can  then  bo  verified  to  be  equivalent  to: 


minimi l*  t  .  E  d„ .y. ,  (10) 

id  J«j'  ij  tJ 

subject  to  the  constraints; 

r  y  *  b.  for  id,  (U) 

r  r  ,  *  s-(  for  iej\  (12) 

id  1 

,  >  0  for  id  and  jeJ*,  and  (13) 

*  d  or  tor  id  and  j«J.  (14) 
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The  problem  (10)- (13)  is  a  standard  transportation  problem  and  banco 
can  be  solved  by  the  primal  transportation  algorithm  (alao  knovn  as  the  MODI 
method  [3]).  We  assume  that  the  reader  is  familiar  with  the  usual  terminology 
that  a  cell  is  an  index  pair  (i,j)  with  row  (source)  i  and  column  (use)  jej; 
a  basis  B  to  the  problem  ( 10) - ( 13 )  is  a  collection  of  (V  +  M)  cells  without 
cycles  (loops  or  stepping-stone  tours)  and  such  that  every  row  is!  and 
column  jeJ  '  has  at  least  one  cell.  A  solution  Cyt j }  1*  basic  if  y^  »  0 
for  (i,j)^  B.  A  basic  solution  is  feasible  if  the  {y^ ^ }  eatisfy  the  con¬ 
straints  (11)- (13).  It  is  well  known  [3}  that  the  MODI  method  yields  a  basic 
optimal  solution  (i.e.,  a  basic  feasible  solution  for  which  cost  Z  la 
minimal)  to  the  problem  (10)-(13). 

DEFINITION  1.  We  define  P  to  be  the  standard  transportation  problem 
(10)-(13)  and  p'  to  be  the  special  transportation  problem  (10)-(14),  We 
define  a  basis  to  be  row-unique,  if  corresponding  to  every  column  jej, 
there  is  an  unique  row  i^  such  that  (i,j)«B  if  and  only  if  i  *  ij. 

By  definition,  B  has  W  +  M  cells.  Since  a  row-unique  basis  has 

s  t 

exactly  one  cell  for  each  column  jej,  it  follows  that  the  M  +  l  column 
has  W  cells;  i.s.,  (.,M+l)sB  for  every  i« 1.  (Such  a  basis  cannot  have 
cycles  since  only  the  M  +  lSt  column  contains  mote  than  00#  coll,) 

Theorem  1  below  establishes  the  connection  between  the  problems  P 

and  P 

THEOREM  1.  There  is  a  ono-to-one  correspondence  between  feasible 
solutions  to  P '  and  row-unique  basic  feasible  solutions  bo  P. 

PROOF.  Consider  any  feasible  solution  fy^J  to  P  \  By  (12)-(14) 
and  from  the  assumption  that  r^  >  0,  It  follows  that  corresponding  to  every 
use  jd  there  is  an  unique  source  ij  such  that  y^  >  0  if  end  ooiy  if 
i  «  i  .  Corresponding  to  this  solution  we  define  B  to  be  the  sot  of  W  y  M 
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cells  {(i.,. i)  Cor  jeJ  U  (i,M+l)  for  id}.  Consequently,  the  solution  {y^} 
is  baric  since  y^.  =  0  for  (i,j)^B»  It  is  a  feasible  solution  for  P 
since  {y.  j  is  feasible  for  P'.  Since  B  is  defined  uniquely,  this 
correspondence  is  unique. 

To  prove  the  converse,  assume  that  we  have  a  row-unique  basic  feasible 

solution  fy  1  to  P.  By  (12)-(13)  and  row-uniqueneaa  it  follows  that  cor- 

ij 

responding  to  every  column  jeJ  there  is  an  unique  row  ij  such  that 

v  =■  r  if  i  =  i  and  zero  otherwise.  Consequently  (14)  is  satisfied 

i.i  J  j 

and  from  (11)- (13 )  it  follows  that  (y^l  is  feasible  to  P'  as  well. 
Furthermore  this  correspondence  is  unique  thus  completing  the  proof. 

By  Theorem  1  and  from  the  fact  that  the  problems  P  and  P  share  a 
common  objective  function  (10)  it  now  follows  that: 

CORO  L  LARI'  1:  There  is  a  one-to-one  correspondence  between  optimal 
solutions  to  P '  and  the  optima  among  the  row-unique  baeic  solutions  to  P. 

A  solution  procedure  to  the  problem  P  now  easily  fellows  somewhat 
along  the  lines  of  the  subtour  elimination  algorithms  for  the  traveling- 
salesman  problem  15,  7|. 

This  algorithm  is  basically  a  branch-aml-bound  procedure  which 
begins  by  partitioning  the  set  of  row-unique  basic  feasible  solutions  and  then 
••alee  lacing  lower  hounds  on  the  costs  of  all  solutions  le  0  subset.  Ihs 
initial  hound  is  found  by  solving  the  standard  transportation  problem  ?. 

U  the  basic  optimal  solution  to  P  is  row-unique  then  we  ere  finished  in 
the  sense  that  we  have  an  optimal  solution  to  P  as  well  (Corollary  1), 
Suppose  .->n  the  contrary  that  the  basic  optimal  solution  to  P  is  not  KSf* 
unique.  Let  us  denote  by  J  '  one  of  the  columns  jcJ  which  has  more  than 
than  one  tell  belonging  to  B  and  let  (i  ,j  )  be  one  such  basic  cell. 
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(Though  any  such  (i'.j')  can  be  chosen,  we  discuss  below  a  heuristic  for 
cnoosvng  a  ’good’  (i  \j  ')  from  the  point  of  view  of  computational  efficiency.) 
Vie  now  branch  into  two  subproblems  (a)  the  subset  in  which  (1  ')  is  a  cell 

in  the  optimum  row-unique  basic-optimal  solution;  and  (b)  the  subset  in  which 
( i  ' ,j ' )  is  not  a  cell  in  the  optimum  solution.  The  two  new  transportation 
problems  corresponding  to  (a)  and  (b)  are  solved  to  determine  the  lower  bounds 
t.ov  all  row-unique  basic  optimal  solutions  in  their  respective  subsets.  If 
the  optimal  solution  corresponding  to  any  one  subset  is  rov-urique  and  the  cost 
of  th* s  solution  is  less  than  or  equal  to  the  lover  bounds  on  all  other  subsets 
then  r.uch  .1  solution  is  optimal.  If  not,  Chen  one  selects  that  subset  having 
Muallosr.  lower  hound  and  branches  again  into  two  subjroblem*.  Eventually 
i>m>  Is  {iHstrmi  of  finding  an  optimum  row-unique  basic  optical  solution  and 
vc  a.'u  quenrly  »n  op  liman  to  V  (by  Corollary  1). 

fM-vyvn  l  e  unwont  e  on  the  above  a  Igor  it  hie  are  now  in  order.  first,  it  is 
c»h*.?hvi-r  that  rhy prorodure  for  branching  on  a  non  row-unique  basts  excludes 
that  b4-:.is  '?r««  t  ho  two  subsets  hue  does  not  exclude  any  tow-unique  basis. 
ih<;  algorithm  converges 'in  a  finite  nussber  of  steps  since  the  total  number 
><?  !.,i  i  tinUe  and  since  at  least  one  basis  is  excluded  at  .every  UttiUwt. 
Second;  th«j  branching  procedure  rueuUe  in  e  partition  of  the  row-unique 
has  it  mlutichs  in  that  subset  and  hence  the  elgorith*  c«a  be  expected 

so  bs  t'l-fleirnt.  ’Jhird,  for  the  subp ruble*  with  constrained  to  be  in 

lh*»  opt  solution,  by  vow-uniqueness  it  follows  that  <1  *)  is  the  only 

*  -5  5  in  wo  Sown  i’.  Consequently,  we  can  drop  column  j#  fro*  further  con* 
'..d»>5 at mm*  no-utv  b  .  to  b.  ,  -  r  *  and  solve  a  seal  Ur  transportation 
pfoblers.  ihu  reduction  in  ,  may  further  simplify  the  proble*  since  the 
r < i  for  which  r^  to  greater  than  the  ae»  vaiua  of  bfc*  cannot 
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possibly  bo  in  the  optimum  soluticn  (such  cells  can  be  eliminated  by  d#» 
fining  d^  .  *  «).  Fourth,  the  optimal  solutions  to  the  aubproblem*  can  be 

.  f f iciently  obtained  by  the  operator  theory  of  parametric  programing 
developed  in  [6,  0)  rather  than  re-solving.  Moreover,  tb*  backtracking  steps 
oi  the  branch  <nd  bound  procedure  may  also  be  done  this  way.  Finally,  a 
non- row-unique  basis  can  have  at  most  W  columns  vhieh  have  more  than  one 
evil  (this  follows  from  the  fact  that  a  basis  has  W  *  K  cells  vtth  at 
least  one  cell  tor  each  of  the  M  columns).  Consequently  the  fraction  of 
ibe  naxtaua  t.  umber  of  columns  which  do  not  satisfy  row- uniqueness  is  V/M. 


Vht.f  the  proposed  algorithm  can  be  expected  to  be  relatively  ®or*  efficient 


tot  problem*  where  the  number  of  uses  greatly  exceed*  Use  number  of  sources. 

We  new  consider  the  question  of  choosing  the  cell  (i  \ j  *)  upon  which 

* 

•o  wake  the  computational  procedure  branch.  Let  u«  denote  by  J  t he  iet  of 

.•.•Iww*.  tbas  have  two  or  more  cell*  of  the  basis  U*C  J) ,  Given  a  colvssa 

w  suggest  branching  on  that  basic  cell  (i  ,j)  for  which  d^  i*  the 


*n*U«*t.  than  along  the  branch  in  which  (i»j)  i*  excluded  from  the  optimal 

«*»luu?n  the  teat  tan  bo  espec.  A  to  increase  approximately  by  ^  *  *d£)"di 

where  ■  ks  the  next  smallest  cost  of  a  basic  cell  le  column  )  and 

i«  ? :sssoun£  shipped  via  the  scaliest  coat  basic  cell  (1, j)«  Couaeqmantly, . 

t  a 

(at  branching,  «•?.  tan  choose  the  column  )  cJ  for  which  dj  la  the  largest 
and  branrh  nti  (»  \  J  *)  -  where  <l  Vj  ’)  ha*,  the  lowest  coat  among  all  baaic 

>: *■  si-  ■  i*;  C'rhiib  |  *. 


ij 


wc  aasf8»4r««c  the  above  results  in  Algorithm  l  for  Mlvlog  the  apacial 

irapspo*  tat  ion  ptew*  **•«. 


AS-GGSSflSK  t.  t  finding  an  optimal  solution  to  tit  special  teaaaporta- 


tlroi  fntMet  fU-ttj. 
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Sot  up  the  problem  P  defined  by  (10)-(i3).  let  denote  the 
problem  P,  «  0  denote  the  set  of  cell-s  constrained  to  be  included 
in  the  optimum  solution  and  «  0  denote  the  of  tells  excluded 
from  the  optimum  solution.  Let  be  the  optimum  solution  to  Pj 

with  basis  dj  and  cost  2^  (this  stay  be  obtained  by  the  print.  (MODI) 
method).  Let  S  =  { 1 }  denote  the  set  of  problems  under  consideration  and 
let  m  -  1  denote  the  total  number  of  problems  generated  so  far. 

Choose  the  problem  for  which  2^  is  the  smallest  for  UeS.  If 
is  row-unique  *o  to  (5).  Otherwise  go  to  (3), 

ta)  Find  t ho  set  of  columns  J  for  which  the  basis  has  two  or 

* 

more  basic  cells  in  that  column.  For  each  column  J«J  find  the 

two  basic  cells  (i.j)  and  (4,j)  for  which  the  unit  coats  ar«  *fce 

smallest  a-wi  the  second  smallest  respectively.  Define  ^  •  (d^-d^ly^ 

and  choose  the  *  *ej*  for  which  a,  ta  the  lar|*»t.  Select  the 

I  vest  :;un  e.*stc  cell  U'.j*'  in  column  J /,  for  branching. 

(h)  nofinv  P_4l  a*  the  problem  obtained  from  by  const rainlnt 

ii  '.'Set'  be  an  addition*!  basic  cell  i.e.,  *  A,.  U  {{i  *,J '')} 

«««  let  .?  j  *  The  prmblco  P  ^  can  be  obtained  itm.  by 

dropping  coJuan  ]*  and  Jo fining  b^  •  to  b#  bj »  -  v y  For  column*  j 

swh  ths*  r  •  b.  »  define  d.  *  *  ®», 

.1  i  »  ) 

(c)  Define  I-  os  the  problem  obtained  from  P,  by  excluding 

■  sfr’Si  K 

i  r»v.  tho  optimal  boils.  Set  ui(f  \j  *}}  H.i.,  d&  #^  #  «  •) 


and  l<rt  f.  .«  M  »L  . 

%H2  «t 


i)  ifc>«otc  the  basic  optimal  solutions  to  Vl  and  (obtained  by 

the  MPhl  sseUwd)  to  bo  Y^j  and  with  bait#  and 


Define  Z^j  «  optimal  cost  to  + 


E  r.d. .  tad 


Z..,  *  optimal  cost  to  P. ,  +  S  ri(, 

"  1  (i.Dca^j  1 11 
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4.  Drop  k  from  Lhe  set  S  and  add  (ro+1)  and  (m+2)  to  S.  Redefine  a 
as  m+2.  Go  to  1.2 j. 

5.  The  optimal  solution  to  the  special  transportation  problem  (1-4)  is  given 
by  Y.  and  0^  with  the  associated  cost  *  optimal  cost  to 

P,  +  I  r  .d.  . .  Stop. 

k  i  1 1 

(i,j)€\  J  J 

We  illustrate  below  the  application  of  Algorithm  1  with  the  same 
example  as  was  solved  in  [4].  Fig.  la  shows  this  problem  with  four  sources 
S1,...,S4,  five  uses  Ul,...,U5,  costs  c^,  capacities  b^  and  demands 
as  shown.  Note  that  since  >  b^  and  r2  >  b^,  the  cells  (1,4)  and 
(2,4.)  cannot  possibly  be  in  the  optimum  solution.  Consequently  ”  cc*. 

Figure  la  -  id  about  here 

At  stop  (1)  of  Algorithm  1  we  set  up  the  transportation  problem  *  P 

as  shown  i.n  Fig.  lb,  by  adding  a  dummy  use  U6  with  demand 
4  3 

v  ~  I  b.  -  I,  r.  -  14-11-3  (eqn.  (9))  and  defining  costs  d  as 

(’  i»l  '  j*\  J  J 

per  equations  (6)  and  (8).  For  the  problem  none  of  the  calls  are 

constrained  to  be  Included  or  excluded  in  the  optimum  solution  so  that 
Oj  =  7 ^  -  0.  The  optimum  solution  to  P1  obtained  by  the  primal  method  (the 
capacities  b .  were  perturbed  slightly  to  prevent  cycling  {3])  is  also  shown 
in  Fig,  la  where  the  circled  cells  denote  the  basic  cells  with  the  amounts 
written  over  the  circles  (y^  =  0  for  non-basic  cells).  The  optimum  value 
for  the  objective  function  can  be  verified  to  be  *  19/3.  We  now  set  tn  *  1 
■and  S  {l}.  In  step  (2)  of  thn  algorithm,  we  find  that  the  basis  of  Fig.  lb 
is  not  row-unique  so  Unit  we  proceed  to  step  (3). 

in  stop  3(a)  wo.  find  .)  ™  {l , 2 }  so  that  *  (2/3  *  1/3)  x  1  *  1/3 

and  0*2  »  (1  -  1/3)  x  2  *  4/3  so  that  j'  *  2  and  (i \j  ')  •  (2,2).  In 


step  3(b)  problem  P^  is  defined  as  problem  P^  with  (2,2)  constrained  to 
be  included  Ln  the  basis  (i.e.,  s  1(2,2) j  and  Yj  *  $)«  Conaequently , 

ue  drop  use  U2  and  change  b0  to  b.,  •  a  4  -  3  *  1.  Since  r^,  and 
are  greater  than  b^,  the  cells  (2,1),  (2,3),  (2,4)  can  not  be  in  the  optimal 
solution  so  that  we  set  d7^  =  =  «  and  obtain  P,,  The  optimal 

solution  to  P  shown  in  Fig.  lc . ,  was  obtained  by  the  primal  method.  The 
optimal  cost  of  the  solution  to  P^  can  be  verified  from  Fig.  lc  to  be  23/3 


so  that  Z_ 


(23/3)  + 


r  ,d 


j  ij 


(23/3)  +  (3  x  1/3)  «  26/3.  Similarly 


is  obtained  from  by  excluding  (2,2)  from  the  optimal  solution  (0.^  * 

=  {(2,2)}).  ioi  equently  d7^  is  sat  equal  to  «  in  Fug.  Id.  The  optimal 
solution  to  P^  is  also  shown  in  this  figure  with  *  23/3.  The  branching 
of  to  P2  and  on  the  basis  of  cell  (2,2)  can  be  seen  in  Fig.  2 

as  well.  We  now  set  S  *  [2,3}  and  m«*3  and  return  to  step  (2), 

Since  Z 3  <  Z^  and  since  is  not  row-unique  we  now  branch  the 

problem  P.}  into  two  subsets.  From  Fig.  id.  J  «*  {l,4},  °  2/3  and 

.1,  =  l  so  that  j'  =  4  and  *  (3,4).  In  step  3(b)  we  define  P^ 

to  be  the  same  as  but  with  (3,4)  constrained  to  be  included  in  the 

optimal  solution  (i.e.,  ~  {(3,4)}  and  *  {(2,2)}).  Consequently  we 

drop  and  change  to  *  3  «  2  »  l.  Since  r j ,  rj>  v-\  are 

greater  than  we  make  d^  «=>  n  d^  =  <*>  to  obtain  Fig.  le.  The  optimal 

solution  to  Is  shov.n  in  Fig.  le  with  cost  25/3  so  that  "  25/3  +  (r^  x d.^)  "  28/3 


Figures  le  -  lh  about  here 


Fig.  1!:  shows  the  problem  P^  obtained  irora  P^  by  constraining  (3,4) 
to  he  excluded  from  the  optimal  solution  (i.e,,  a  {(2,2),  (3,4)}),  We 
mark  *  »  and  obtain  the  row-unique  basic  optimal  solution  of  Fig.  If.  with 

Z,.  ■*  9,  In  .step  (4),  ii  becomes  {2,4,5}  and  m  «*  5. 


TfPCf'C JTF^> 
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Wc  now  return  to  step  (2)  of  the  algorithm  to  find  that  Z2  ia  the 

smallest  among  the  problems  in  S  so  that  we  branch  P  to  P,  and  P 

It  1 

on  the  basis  of  cell  (3,1)  as  shown  in  Figs.  Ig  and  lh.  Now  S  becomes 
[4, 5, 6, 7}  so  that  is  the  smallest  cost.  Since  ia  row-unique, 

the  optimal  solution  to  the  special  transportation  problem  is  given  by  Fig.  if 
with  cost  Z5  =  9.  This  optimal  solution  assigns  the  uses  Ul,  U2,  U3,  U4,  US 
to  sources  S3,  SI,  S2,  S4  and  SI  respectively,  the  same  solution  as  in  {4], 
Figure  2  shows  the  branch-and-bound  tree  at  the  end  of  the  computation. 

Figure  2  about  here 


From  a  computational  point  of  view,  it  is  not  necessary  to  score  the 

problems  F^  for  keS.  It  is  enough  if  we  store  the  sets  Q,  and  Y  for 

keS.  To  construct  P^  from  the  original  problem  P,  wo  first  set  *  ® 

for(i,j)eY.  Next,  for  every  (i  ,j)  e  Q.  we  drop  column  \  and  modify  b. 

K  1 

to  bi  *  rj’  Finall>'  wp  eliminate  those  cells  (i,j)  with  jej  for  which 
tv  >  b .  (by  defining  d^.  “  «). 

It  is  interesting  to  compare  our  algorithm  to  the  implicit  enumeration 
approach  in  [4],  The  latter  starts  out  with  the  solution  of  the  total  cost  Z' 
obtained  when  each  use  j  is  assigned  to  that  source  i  with  the  least  coat 
;v  . ,  The  feasibility  condition  (14)  is  satisfied  at  every  step  but  not  (11). 
On  the  other  hand,  our  procedure  starts  with  the  least  cost  optimal  solution 
to  P  of  cost  Z  "  and  maintains  the  feasibility  condition  (11)  but  not  (14). 
Denoting  by  Z  the  optimal  cost  of  the  special  transportation  problem,  the 
relative  efficiency  ot:  the  algorithms  will  vary  across  problems  depending  on 

i  i  f  it 

whether  Z  or  7,  is  closer  to  Z  .  As  mentioned  earlier,  for  problems  with 


iar 


go  .M/ u  the  i  hi  c  ■*.*;  i  b  I  li  ly  .if  (14)  is  relatively  small  so  that  our  al¬ 


gorithm  is  better  suited  lor  such  problems .  On  the  other  hand  for  problems 
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;:cr  which  the  ratio  M/W  is  small,  the  all  zero-one  algorithm  of  {&}  can 
he  expected  to  be  more  efficient. 

Finally  it.  should  be  pointed  out  that  although  this  algoritho  has 
been  developed  using  the  primal  method  as  a  subroutine  for  solving  transporta¬ 
tion  problems,  other  methods  such  as  the  primal-dual  methods  could  also  be 
used . 


2.  EaTEF  10N„  AND  APPLICATIONS 

We  first  formulate  an  extension  of  the  special  transportation  problem 

where  th<  capacities  h.  can  to  increased  by  unit  cost  Demoting  by  the 

additional  -apaclty  of  source  i,  equations  (10)  and  (11)  arc  «oci tied  to 


become  (.Ibl-  'lb)  bei.au*: 

A  n  ~  Ur«,  +  Z  d  y  ]  and  (15) 

id  ‘  1  JCJ  '  lj  41 

r  y.  «  b  +  u.  for  Ut  (16) 

i  c.i '  J 

Lot  ns  denote  by  h,  the  maximum  additional  capacity  that  ran  be  added  to 
sour co  l  i  if  time  Is  no  -oen  constraint,  h.  can  be  sec  equal  to  a  very 
large  number ) .  As  a -further  genet*.  .‘sac  ion  let  ^  denote  the  unit  coat 
oi  not  utilizing  the  capacity  oi  se  .veo  t  (if  this  Involves  *  unit  .eving 
then  .p.  would  hr  .negative)  ant  lot  q^  (  j.  ®)  denote  Che  wintenaa  Mtllita- 
turn  level  for  source  l.  Defining  q^  ■  b^  -  q^  the  following  relations 
.hold  l "t‘  tin*  slack  into  M  +  I: 


i ,M+1  "  Pt 

lor 

id,  and 

u  n 

i ,M+1  - 

for 

Wl. 

08) 

The  additional  capacities  can  b.u  thought  oi  as  a  aurplua  use  (H  4  2), 


We  now  define 
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J  "  =  J  '  U  {  (M+2  )  ! ,  a  ml 

y.  a  h.-u.  ter  t. el . 

'i,M+2  1  1 

Furthermore,  since  0  <  •_£  h^.  we  have 

0  —  -  **l  toriCl' 

The  objective  function  (15)  becomes 


Z  =  s0  + 


id  jsj 


T'  <1.  .y.  . 

-  tj  tj 


when 


ZA  =  ?  &;b.  and 

id 


J1,H«  ’  *»l  l0‘'  UU 


The  constraints  (16)  become 


f  y.  ,  *  b.  4-  h.  for  icl . 
Hi*  Lj  1  1 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 


Hit’  const raints  (12)  nold  as  usual  for  jcJ.  But  for  tho  dummy  usee  (M+l),  (12) 
Should  be  modified  to 


.  “  'i.  t-H  .  . 
id  id 


v  +  u4 >  ”  t  r , 
jeJ  •' 


(b.  +  h.)  -  ~  r 


J 


(26) 


its*  that  from  (20)  wo  obtain 

*2  y.  ,,.«•<*  £  y.  u ,  ^  *  \vj  * 

U*t.  i,N  1  id  i,H “  id  1  i  jeJ 

The  etMUJt i aint  (26?  is  nor  a  regular  transportation  constraint  since  it 

involves  variable**  from  two  columns.  To  bring  it  to  standard  transportation 

term  we  deft  tic 


"SL  •  ,rt  6,h*i  an0 

l  51 


Vv.'+1  ,M+2  “  n2  ‘  Ul  yi.M+2 


(27) 

(28) 


where  Nj  and  N.,  are  largo  positive  numbers  so  that  y^+^  and 
v(,^.  are  noimegativv.  Consequently  (26)  becomes 
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Figures  3  summarises  the  capacitated  (or  upper  bounded)  trunnpi'ttaciou 
formulation  of  thia  problem.  The  special  transportation  problem  has  the 
additional  conatraint  (14)  that  each  use  jeJ  has  to  bo  supplied  in-  only 
one  (poaaibly  different)  source  id. 

An  algorithm  for  this  generalized  problem  should  bo  obvi  m.s,  Wo 
can  utilize  the  tame  branch  and  bound  procedure  of  Section  2  with  the 
capacitated  transportation  formulation  ot  Figure  3,  However,  ehe  implicit 
enumeration  algorithm  of  [4]  is  not  capable  of  such  an  easy  extension  (ol- 
though  DeMalo  and  Roveda  (4]  in  their  concluding  discussion  suggest  the 
problem  generalisation  considered  here). 

Though  the  special  transportation  model  concerns  itself  with  sources 
and  uses  typically  considered  as  warehouses  and  markets,  wo  wish  to  point  ou 
that  it  offers  en  Important  generalization  to  assignment  models.  (Km  other 
interesting  and  important  assignment  problem  generalization.-  ,n-c  the  paper  ( 
by  Charnes,  Cooper,  Nichaus  and  Stedry.)  Consider,  for  instance,  assigning 
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jobs  Co  machines  in  a  case  when  it  may  be  prohibitive  Co  do  the  same  job 
on  mere  than  one  machine  (perhaps  because  of  set-up  cost  considerations). 
Denoting  by  r^  the  time  required  to  perform  the  job  J,  b^  the  time 
available  on  machine  i  and  the  cost  of  performing  job  j  on  machine  i, 

we  obtain  the  special  transportation  problem  (l)-(4).  Similarly,  this  model 
can  also  be  utilized  in  assigning  workers  to  supervisors  (or  students  to 
advisors)  where  r^  is  the  time  needed  to  supervise  the  j-th  worker.  These 
applications  suggest  a  further  extension  of  problou  (l)-(4)  where  (2)  is 
replaceu  by 


E  r. ,  x, .  <  b.  for  iel  , 
J.J  li'  * 


(34) 


i.e.,  r^  ia  not  necessarily  constant  for  all  id;  in  other  words  job  j 
might  be  done  with  differing  efficiencies  by  each  of  the  machines.  Tha 
branch  and  bound  procedure  of  Section  2  would  then  have  to  be  applied  to  a 
generalized  transportation  problem  [1,6]  with  column  demands  equal  to  unity. 


4,  CONCLUSIONS 

In  this  paper  we  have  considered  a  special  class  of  transportation 
problems  of  assigning  uses  to  source*  and  provided  a  branch  and  bound  solution 
procedure  with  the  standard  transportation  problem  as  a  subroutine.  Compared 
to  the  implicit  enumeration  approach  in  [4]  this  algorithm  appears  to  be  com¬ 
putationally  more  efficient  particularly  for  problems  where  the  number  of  uses 
greatly  exceeds  the  number  of  sources. 
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Note  l 
Note  2 


Fig.  2 

Branch  and  Bound  Tree  Diagram  at  Optimum 

R  ■  Row-unique  optimal  basis;  NR  ■  Non  row-unique  optical  basis 
The  label  (2,2)  connecting  and  Pj  indicates  that  is  obtained 
from  by  constraining  its  optimum  solution  to  include  the  cell 
(2,2).  Similarly  is  obtained  from  Pj  by  excluding  the  cell  (3,4) 
from  its  optimal  solution  (denoted  by  (3,4)). 


Uses 


Sources 


Dummy 

Source 


Slack  Surpl  ■« 
use  use 


"i  “a  •  •  •  U ,  •  •  •  Vl  V.'  Capacities 


Demands  rl  •  '*2  *  ’ ' *  ri 


i 

Transportation  Format  lor  the  Generalised  Problem 

Note;  l.  In  each  cell  the  number  at  the  center  denotaa  the  unit  c«at  d^» 
the  number  at  the  opper  righchand  corner  denote*  an  upper  bound 
for  the  cell  ii£  this  is  blank  this  implies  that  there  is  no 
copper  bound), 

2,  (i,  K^,  denote  very  large  positive  nuabers. 


